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Abstract. Gravitational lensing time delays depend upon the Hubble constant and the density distribution of the lensing galax- 
ies. This allows one to either model the lens and estimate the Hubble constant, or to use a prior on the Hubble constant from 
other studies and investigate what the preferred density distribution is. Some studies have required compact dark matter halos 
(constant M/L ratio) in order to reconcile gravitational lenses with the HST/WMAP value of the Hubble constant (72 ± 8 km s _I 
Mpc~' and 72 + 5 kms~' Mpc -1 , respectively). This is in direct contradiction with X-ray, stellar dynamical, and weak lensing 
studies, which all point towards extended halos and isothermal density profiles. In this work, we examine an up-to-date sample 
of 13 lensing galaxies resulting in a data set consisting of 21 time delays. We select systems in which there is a single primary 
lensing galaxy (e.g. excluding systems undergoing mergers). Analysis is performed using analytic models based upon a power- 
law density profile (p oc r^) of which the isothermal profile is a special case (tj = 2). This yields a value of t] = 2.1 1 ±0.12 (3cr) 
for the mean profile when modeling with a prior on the Hubble constant, which is only consistent with isothermality within 3 
cr. Note that this is a formal error from our calculations, and does not include the impact of sample selection or simplifications 
in the lens modeling. We conclude that time delays are a useful probe of density profiles, in particular as a function of the 
environment in which the lens resides, when combined with a prior on the Hubble constant. 

Key words. Gravitational lensing - Cosmological parameters, dark matter, halos 



1. Introduction 

In 1964, Refsdal first suggested that time delays between gravi- 
tationally lensed images could be used to determine the Hubble 
constant Hq. The appeal of this method was the fact that it did 
not depend upon the local distance scale and was a relatively 
straightforward single-step process. In addition, the time delay 
theory upon which it is based was derived from the established 
theory of general relativity. Only in recent years, however, with 
the availability of high resolution radio and optical imaging, 
has the method begun to yield results. 

When focusing on individual lensing galaxies, a large num- 
ber of lens modeling studies (e.g. Impey et al. 1998; Williams 
& Saha 2000; Winn et al. 2002) have obtained values of H that 
have been lower than that from the Hubble Space Telescope 
(HST) Key Project (Freedman et al. 2001) which yields a 
value of Ho = 72 + 8 kms^Mpc -1 , or indeed the Wilkinson 
Microwave Anisotrophy Probe (WMAP) which yields a very 
similar Hq - 12 ± 5 kms^'Mpc -1 , although this assumes a 
flat universe (Spergel et al. 2003). In particular, a few studies 
have inferred that if one assumes that lenses have an isother- 
mal mass distribution (flat rotation curves) the derived value 
is Hq ~ 50kms _1 Mpc _1 , whereas if they assume a constant 
mass-to-light ratio (M/L), which results in a falling rotation 
curve, then Hq ~ 70kms~ 1 Mpc~ 1 , which is in agreement with 



the HST/WMAP finding (Kochanek 2002; others e.g. Impey et 
al. 1998; Rusin et al. 2003). This constant M/L ratio implies a 
compact dark matter halo. 

However, other studies such as weak lensing (e.g. Guzik & 
Seljak 2002), stellar dynamics (Rix et al. 1997; Romanowsky 
& Kochanek 1999; Gerhard et al. 2001; Treu & Koopmans 
2002), and X-ray galaxy measurements (e.g. Fabbiano 1999; 
Lowenstein & White 1999), are consistent with near isothermal 
density profiles (i.e. extended dark matter halos). Note however 
that in the strong lensing regime we are probing the inner few 
kpc of galaxies (i.e. within the Einstein radius), whereas other 
studies, such as galaxy-galaxy lensing, are sensitive to more 
extended scales. 

The aim of this paper is to study whether time delays 
from multiple image lens systems statistically point towards 
extended dark matter halos which are inferred from other mea- 
surements, when an external prior is adopted for Hq. We per- 
form the study without the complexity of detailed lens mod- 
eling in an effort to provide insight with a simple model. The 
outline of this paper is as follows. In §2 we introduce two ana- 
lytic models for the gravitational time delay based upon an SIS 
potential and general power law potential. §3 will discuss the 
effect of ellipticity and shear on the time delays. In §4 we dis- 
cuss the lens sample used as the data set in the investigation, 
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and in §5 we explain the analysis of this data set and present 
results. §6 draws some conclusions based upon the findings. 

2. Analytic Models for the Gravitational Time Delay 

The light-travel time for an image of a gravitational lens sys- 
tem, is given by (see e.g. Narayan & Bartelmann 1997) 



M 



l+z. 



D,D S 



D, 



l -\o-p\ 2 -m 



(i) 



where zi is the redshift of the lensing galaxy, D, and D s are the 
angular diameter distances to the lens and source respectively, 
and D[ s is the angular diameter distance between the lens and 
source. The two dimensional vectors 6 and /?, are the (angu- 
lar) image and source positions respectively and the effective 
potential at the locations of images is i[>{0). Note that the light- 
travel time has a geometrical component and a component due 
to gravitational time dilation (Shapiro delay) which depends 
on the gravitational potential. Since angular diameter distance 
scales as c/Ho, where Ho is the Hubble parameter at the present 
epoch, it follows that t oc \/Hq. The time delay between two 
images, say A and B, is the difference in their light-travel time 
i.e. At — tA - ts- In order to calculate this quantity we must 
adopt a specific effective potential. 




Fig. 1. Diagram of a time delay lens with two images, A and B. 
The lensing galaxy is at the centre with the images at radii \0\ 
which we denote by Ra and Rb- The images bound an annulus 
of width AR - R A - Rb, with average radius (R) = (R A +Rb)/2 
(as in Kochanek & Schechter 2004). 



2.1. Time Delay from an SIS Potential 

The singular isothermal sphere (SIS) is a simple yet surpris- 
ingly realistic starting point for modeling lenses. Its lensing 
potential has the form: 



m = b\ 



(2) 



where 



b = 4n- 



D s c 2 



is a deflection scale determined by the geometry. For a < b, 
the SIS produces two collinear images at radii Ra = \P\ + b and 
Rb — b - |/3| on opposite sides of the lens galaxy, as in Fig.^ 
The time delay between these two images can be shown to be 
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which depends upon the annulus bounded by the two image 
positions, as in Fig.^ and as derived in Witt et al. (2000). 

2.2. Time Delay from a General Power Law Potential 

A more general approach is to use an effective potential that 
allows for different slopes (77) for the radial density profile 



3-D 



of which the SIS is a special case (77 = 2), with more cen- 
trally concentrated mass distributions having higher values of 



77 (tending to a point mass with 77 = 3). The deflection scale 
b is given by (e.g. the contribution of Kochanek in Kochanek, 
Schneider & Wambsganss 2004) 
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and the convergence is 
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Previous studies have found that lens galaxies are well de- 
scribed by such power-law profiles, in which p oc r~ n (e.g. Witt 
et al. 2000; Rusin et al. 2003). To obtain a time delay from 
the power-law lens potential, we use the expression detailed in 
Kockanek (2002): 
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is the mean surface mass density in the annulus bounded by 
images A and B, in units of the critical surface mass density. 
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In the analysis that follows in §5, we use £] M = 0.3, f2 A 
= 0.7, Q. K = 0, which enter through the angular diameter dis- 
tances. 

3. The effect of Galaxy Ellipticity and External 
Shear 

In the preceding section we introduced a spherically symmet- 
ric effective potential which will be used to model the lensing 
galaxy of all systems in the data set, including four-image sys- 
tems. Our goal is to obtain an estimate of the average slope of 
the density profiles using a simple approach, comparing sys- 
tems in a homogenous manner and without detailed modeling 
of ellipticity or external shear. We will now, however, briefly 
discuss the effect these have on the corresponding time delays 
using Fermat surfaces for illustration. 

Fig.|2]shows four Fermat surfaces, as derived from equation 
0, each of which displays a different image configuration de- 
pendent on the effective potential and parameters used. For all 
four plots the source position, velocity dispersion, lens/source 
redshifts and Hubble constant were kept fixed. We consider the 
changes in time delays of the images A, B and C, present at the 
stationary points of the surface (note that there may be other 
stationary points). We investigated values of ellipticity up to 
0.2, this being a moderate ellipticity for a lens galaxy. 

Fig-Et begins with a standard SIS potential, with a time de- 
lay of AtAB = 51.2 days. An ellipticity is then introduced into 
the effective potential (an SIE with e = 0.1), as seen in (b), giv- 
ing a time delay of AtAB = 57.1 days. Finally, in plots (c) and 
(d), we consider the impact of changing the slope of the effec- 
tive potential. Keeping e = 0.2, the slope of the density profile 
is allowed to change from 77 = 2.0 to 77 — 2.3 (corresponding to 
a slope of 1.0 to 0.7 in the effective potential). The time delay 
is AtAB = 67.9 days and Af^c = 48.8 days for the 77 = 2.0 case, 
and AtAB = 67.2 days and Af^c = 46.8 days, for the 77 = 2.3 
case. In the limit 77 tends to 3.0, AtAB = 66.7 days and AtAc = 
44.5 days . 

From this we highlight two key features: firstly, the time 
delays between images A and B change by less than 10% when 
comparing the SIS to the SIE with a small ellipticity. Secondly, 
when we consider the change in the time delays as 77 goes from 
2.0 to 2.3, keeping e constant, we observe very little change 
at all. We pick these particular values for illustrative purposes 
only. Indeed, moving beyond 2.3 to higher values of 77 had a 
similarly small effect. 

As noted by Witt et al. (2000) for a generalised isother- 
mal lens (allowing for general angular structure irrespective of 
ellipticity) the Hubble constant is simply related to the image 
positions and time delay along with the source and lens red- 
shifts, not requiring detailed models. When the profile deviates 
from isothermal, and ellipticity is included, detailed modeling 
cannot be avoided. However, for images which are not close 
pairs, they derived scaling relationships which show that for 
moderate ellipticities, the time delay is well approximated by 
the isothermal delay modified by a factor which depends on 
the slope of the density profile or effective potential. The scal- 
ing relations themselves follow the form (77- 1) for images at 
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SIS SIE (e = 0.1) 
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Power Law + Ellip. (e = 0.2, ri = 2.0) Power Law + Ellip. (e = 0.2, n = 2.3) 
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Fig. 2. A plot of four Fermat time delay surfaces for varying 
effective potentials. The images are to be found at the station- 
ary points of the surface, and are displayed as A, B, and C. 
The corresponding time delays in units of days are as follows; 
(a) Atsis =51.2; (b) Ar s/£ . £=0 .i =57.1; (c) Af, ;=2 . , e =o.2 = 67.9 
(AB), Af, ;=2 . , £= o.2 = 48.8 (AC); (d) A/,, 2 . 3j , .2 = 67.2 (AB), 
A^ =2 .3, e =o.2 = 46.8 (AC) 



~ 180°, and (T]-l)/(3-rj) for those at ~ 90°, and are multiplied 
with the isothermal time delay Af s , s , of equation ©. 

To estimate the effect of these scalings we generated image 
pairs (at both ~ 90° and ~ 180°) from Fermat surfaces derived 
from the parameters e = 0.2 and 77 = 2.3 (see Fig.|2jl). Fixing 
//o = 72 km s _I Mpc~', we recovered 77 from the time delay 
expression in §2.2 (eqn. |8j for both sets of image pairs. We 
found that when using the 180° scaling we recovered a value 
of 77 = 2.32 for image pair AC, and 77 = 2.46 for AB. Using the 
~ 90° scaling yielded 77 = 2.18 for AB. In the limit 77 tends to 
3.0, the maximum error in recovered 77 is less than 20%; note 
that using equation l|8j as we do below, rather than taking into 
account 90° scaling is actually more accurate in recovering the 
slope. 

When considering the effect of the ~ 90° scaling on the real 
image data of Table 1, we found that in general the scaling gave 
values of 77 both higher and lower than those derived from the 
~ 180° scaling. Moreover, the difference in the derived values 
for 77 from both scalings was less than 7%. 

From the above, we see that while using a spherically sym- 
metric effective potential or mass profile in the analysis of four 
image systems is certainly to be considered an approximation, 
it would appear not to be grossly inaccurate when considering 
the variation of derived parameters. 

External shear also impacts on the analysis, as discussed 
by Witt et al. (2000); they show that if there is no knowledge 
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of the angle between the shear and the lensing galaxy's major 
axis, then a 10% shear would lead to a 10% uncertainty in time 
delay. However, strong external shear would lead to more ex- 
treme uncertainties, which will lead to the exclusion of certain 
systems from the working data sample. 

4. Data: The Lensing Sample 

To date, there are 13 gravitational lensing systems with well 
measured time delays. However, a number of these have more 
than two bright images. In total, there are 21 time delay mea- 
surements, and it is these that constitute the initial data set for 
our investigation (see Table 1). The primary source of errors 
in the data originate from the time delay; most of the errors in 
image positions are on the order of a few percent, and as such 
are negligible in comparison. 

In the analysis that follows in §5, there are three lens 
systems that have been excluded from the working data set. 
These are RXJ091 1+0551, Q0957+561 and B1608+656. For 
RXJ091 1+0551, the main lensing galaxy has a smaller, less 
massive galaxy inside the Einstein ring, and is also embedded 
in a large X-ray cluster (Morgan et al. 2001). Q0957+561 is a 
special case as the lensing galaxy is a BCG (brightest cluster 
galaxy). In both cases, there are strong external perturbations 
that would require a much more complex effective potential 
to successfully attempt to model the time delay, and as such 
would not be well modeled by the expressions in §2. Weaker 
external shear results in far less uncertainty in the modeled time 
delays, as discussed in §3. An additional system that has no- 
table complexity is B1608+656. Although simple lensing po- 
tentials appear to produce an acceptable fit (Fassnacht et al. 
2002) we exclude it from the sample considered since the lens- 
ing galaxy appears to be undergoing a merger with a smaller 
galaxy. In addition to the above, the time delays PG1 1 15+080 
(B-C) and B 1422+23 l(A-C) have been excluded since there is 
redundancy with the other delays quoted for the particular sys- 
tem. While each lens system obviously has only one true value 
for the slope of its density profile, each delay that we use from 
that system is an independent probe of that slope, with corre- 
sponding error. 

Finally, we are left with 10 systems and a total of 14 time 
delays for use in the analysis. Note that the systems in the 
lens sample of the Kochanek (2002) semi-analytic study are 
PG1 115+080, SBS1520-530, B1600+434, and HE2149-2745, 
and this sub-sample will also be examined in the next section. 

5. Analysis and Results 

The analysis of the data sample was performed with equation 
of §2. For this simple analysis we exclude the effects of 
ellipticity and shear which allows us to use the time delay ex- 
pression As mentioned in §3, inclusion of either would re- 
quire detailed modeling of each system individually defeating 
the purpose of a simplified approach. For a particular lens sys- 
tem, we use a single component lens model that accounts for 
the sum of the dark and luminous components, and determine 
the time delay for particular values of Hq and tj given the image 
positions along with source and lens redshifts. These values can 
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Fig. 3. The A^ 2 limits for Hq vs. 77 for the time delay system 
PG1 1 15+080 (A-B), including a prior on Hq from the HST Key 
Project. Note the clear degeneracy that exists between Ho and 
77, where more centrally concentrated mass distributions (larger 
77) will result in time delays that produce a higher value for Hq 



then be compared to the observed delay using the^ 2 statistic as 
follows: 



2 (AtuMod - ktuobs) 

x = — 



cr. 



(10) 



Obx 



where Atij y Mod is the modeled time delay from equation 
Atijfibs is the observed time delay for a particular system, and 
cr 2 obs is the error in that observed time delay. This process al- 
lows us to build up a grid of x 1 values in a 2-D plane of Hq 
and 77. As an additional constraint, we add in a gaussian prior 
to expression (II Oi based on the HST Key Project value Hq — 
72 + 8 kms- 1 Mpc" 1 : 



2 = (72 -H ) 2 

X prior 



8 2 



(11) 



One can then relate these x 1 values to confidence levels 
based upon the number of degrees of freedom. Note that when 
dealing with a single time delay there is 1 data point (Ndp = 
1) and 2 free parameters (N/ p = 2), and hence the number of 
degrees of freedom are -1 (since Ndof=Ndp-Nf p ). This means 
that the problem is under-constrained in this instance. 

Despite the inclusion of the HST Hq prior, there exists a 
notable degeneracy in the Hq - tj plane and we highlight this in 
Fig.|3]for the time delay PG1 1 15+080 (A-B) . The direction of 
the degeneracy is such that the time delay data are consistent 
with a range of more (less) centrally concentrated mass distri- 
butions and higher (lower) values for Hq. We include the best 
fit 77 and corresponding l<x error for all the systems included in 
our sample (see Table 2). 
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TABLE 1 

Observational Data for Time Delay Lenses 



Lens/Components 




ll 


Zj 


Ar,j (days) 


Ri (arcsec) 


i?,(arcsec) 


References 


PG 1115+080 (A-B) 





31 


1.72 


11.7 ± 1.2 


1.147 


0.950 


1 


PG 1115+080 (B-C) 





31 


1.72 


25.0 ± 1.6 


0.950 


1.397 


1 


PG 1115+080 (A-C) 





31 


1.72 


13.3 ± 1.0 


1.147 


1.397 


1 


HE 2149-2745 





50 


2.03 


103 ± 12 


0.344 


1.354 


2 


HE 1104-1805 





73 


2.32 


161 ±7 


1.099 


2.095 


3 


SBS 1520+530 





72 


1.86 


130 ± 3 


0.385 


1.207 


4 


Q0957+561 





36 


1.41 


417 ± 3 


1.036 


5.220 


5 


B 1422+231 (A-B) 





34 


3.62 


1.5 ± 1.4 


1.014 


0.961 


6 


B1422+231 (B-C) 





34 


3.62 


8.2 ± 2.0 


0.961 


1.056 


6 


B 1422+231 (A-C) 





34 


3.62 


7.6 ±2.5 


1.014 


1.056 


6 


B 1600+434 





42 


1.59 


51 ±4 


0.250 


1.140 


7 


RXJ 0911+0551 





77 


2.80 


146 ±4 


0.866 


1.327 


8 


B02 18+357 





68 


0.96 


10.5 ± 0.2 


0.100 


0.240 


9 


PKS 1830-211 





89 


2.51 


26 ±4 


0.320 


0.670 


10 


B 1608+656 (B-A) 





63 


1 .39 


76 ± 10 


1.091 


1.085 


1 1 


B 1608+656 (B-C) 





63 


1.39 


30 ±7 


1.091 


0.682 


11 


B 1608+656 (B-D) 





63 


1 .39 


36 ±7 


1.091 


1.451 


1 1 


HE 0435-1223 (A-B) 





45 


1 .69 


8.0 ± 0.8 


1.298 


1.168 


12 


HE 0435-1223 (A-C) 





45 


1.69 


2.1 ±0.8 


1.298 


1.065 


12 


HE 0435-1223 (A-D) 





45 


1 .69 


14.4 ± 0.8 


1.298 


1.302 


12 


FBQ 0951+2635 





24 


1.25 


16.0 ± 2.0 


0.221 


0.879 


13 



References. - (1) Barkana 1997; (2) Burud et al. 2002a; (3) Ofek & Maoz 2003; (4) Burad et al. 2002b: 
(5) Kundic et al. 1997; (6) Patnaik & Narasimha 2001; (7) Burud et al. 2000; (8) Hjorth et al. 2002; 
(9) Biggs et al. 1999; (10) Lovell et al. 1998; (1 1) Fassnacht et al. 2002; (12) Kochanek et al. 2006; 
(13) Jakobsson et al. 2005 



TABLE 2 

Best Fit tj Parameter for Data Sample 



Lens/Components 


Best Fit i] (inc. prior) 


lcr error 


PG 11 15+080 (A-B) 


2.25 


+0.30/-0.27 


PG 1 115+080 (A-C) 


1.92 


±0.18 


HE 2149-2745 


2.56 


+0.44/-0.41 


HE 1104-1805 


1.77 


±0.13 


SBS 1520+530 


2.40 


+0.30/-0.25 


B 1422+231 (A-B) 


1.63 


+0.92/-0.63 


B 1422+231 (B-C) 


2.89 


+0.11/-0.34 


B 1600+434 


2.21 


+0.32/-0.24 


B02 18+357 


2.60 


±0.30 


PKS 1830-211 


1.88 


+0.40/-0.28 


HE 0435-1223 (A-B) 


1.70 


±0.12 


HE 0435-1223 (A-C) 


1.67 


±0.15 


HE 0435-1223 (A-D) 




-/-* 


FBQ 0951+2635 


2.27 


+0.43/-0.32 



* Contours failed to close in the range 1< t\ < 3, yielding no error bounds. 



The result of combining the constraints from all 14 time de- 
lays, including the HST Ho prior, is shown in Fig.0] Since we 
have a total of 14 data points and 2 free parameters, there are 
12 degrees of freedom . The overall reduced^- 2 for the data dis- 
played in Fig.@]is x^ ed - 49.1. Such a high value for a combi- 
nation of all the data reflects the simplified lens model, and the 
lack of a common density profile slope (i.e. that galaxies have 
a range of slopes for their profiles). The extent of the contours 
is greatly reduced over that of Fig.|5J thus providing bounds 
on 77. The best fit density profile slope is r\ — 2. 1 1+0.12 (3<x). 
If instead we exclude redundant time delays from our sample 
based on the largest fractional error rather than absolute error, 
we obtain the same result. This is also the case if we include ev- 
ery time delay, not excluding any which could be derived from 
the others. The selection of our sample described in Section 4 
above means that these results are relevant to lens systems in 
which there is a single primary lens - in other words to galaxies 



which are not undergoing mergers, or in very dense environ- 
ments. Probing the density profiles of galaxies in complex en- 
vironments requires detailed modeling to account for external 
perturbers. Note also that our error bars are formal errors from 
the calculations, and do not address the systematics of sample 
selection or simplifications in modeling. 

When the data are separated into time delays resulting from 
quad image systems and doubles, then from the 7 time de- 
lays coming from each we obtain r\ = 1.79+0.15 (3cr) and 77 
= 2.23+0.2 (3cr) respectively. 

When considering just the four systems PG1 115+080, 
SBS1520+530, B1600+434, and HE2149-2745, before appli- 
cation of the prior on the Hubble constant, we find a steeper 
best fit slope than for the remainder of the sample, in keeping 
with Kochanek (2002). Kochanek et al. (2006) note that this 
may well be due to the environment of the lensing galaxies 
(e.g. halo stripping of group satellite galaxies). Indeed, other 
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2.4 




1 S i I 

60 65 70 75 80 

H„(kms"' Mpc"') 

Fig. 4. Enlarged plot of the combination of confidence limit 
(such as Fig.|3 for all 14 time delays, all of which include 
the HST Hq value as a prior. Plotted are the 68%, 95% 
and 99% likelihood contours. The centroid value corresponds 
to a density profile slope of r\ — 2.11+0.12 (3cr) at Ho = 
71kms~'Mpc~ I . Inset: same plot but on the scale of Fig.[3] 
While the degeneracy is still apparent, the extent of the con- 
tours is greatly reduced over that of Fig. [3] thus providing 
bounds on 77. 



authors have noted that lens galaxy environments could poten- 
tially have a great impact on halo profiles (Dalai & Watson 
2004; Keeton & Zabludoff 2004). This also suggests that the 
profile slopes of the other systems we have excluded from our 
analysis due to the primary lens being in a very complex en- 
vironment would be likely to increase the overall dispersion. 
We also note that if the 4 systems noted above are in addition 
excluded from the analysis then we obtain 77 = 2.06+0.14 (3<x) 
including the prior on the Hubble constant. 

6. Discussion and Conclusions 

Gravitational lens time delays directly depend upon the Hubble 
constant and the density distribution of the lensing galaxies. 
This allows us to either model the lens and infer the Hubble 
constant, or to marginalize over the Hubble constant derived 
from other studies and determine the preferred density distribu- 
tion. Our goal has been to estimate the typical slope of a mas- 
sive galaxy's density profile, using simple models as in Witt et 
al. (2000). 

Although some studies have indicated that it does not seem 
possible to reconcile an isothermal model for galaxy halos with 
the HST/WMAP value for the Hubble constant, rather requir- 
ing a compact halo (Impey et al. 1998; Kochanek 2002), our 
analysis indicates that, on average, this is not the case. For 



the sample of relatively isolated galaxies which we consider, 
as defined in Section 4, the mean value was found to be 77 = 
2. 1 1 ±0.12 (3(f), when incorporating a prior on the Hubble con- 
stant from the HST Key Project. The quoted error arises for- 
mally from the calculation outlined in Section 5 and does not 
include any bias due to sample selection, or due to the simplifi- 
cations of lens modeling - for example, not fitting for ellipticity 
or external shear. Our finding is in agreement with a study of 22 
lens galaxies using self-similar mass models, which resulted in 
77 = 2.07 ± 0.13 (lcr) (Rusin et al. 2003). The high reduced* 2 
of our best fit to the data is a combination of the large amount 
of scatter in the density profile slopes of the sample, and over- 
simplistic models for the lenses. The steeper slope obtained for 
lenses resulting in double image configurations is perhaps sur- 
prising in that multiple imaging is more effective for steeper 
profiles; however, the cross-section for quads scales roughly as 
e 2 , so this also a competing factor. Previous studies (e.g. Witt et 
al. 2000; Kochanek 2002) have shown that two-image systems 
are relatively insensitive to quadmpole structure, and above we 
briefly considered how deviations from SIS impact on the time 
delay surface and estimates of the slope of the density profile. 

We also note that the compact halo that tends to 77 = 3 
is not typical. However, when we consider the four systems 
PG1 115+080, SBS1520-530, B1600+434, and HE2149-2745, 
there is a tendency towards a steeper than isothermal density 
profile, although in this case we are dealing with small num- 
ber statistics with only four systems to constrain the 77 param- 
eter. Kochanek et al. (2006) have recently pointed out that at 
least three of these four systems might have steeper density 
profiles as a consequence of environment - the distinction be- 
ing whether at one extreme the galaxy is a satellite galaxy (and 
partially stripped of its dark matter halo) or at the other extreme 
is a central galaxy in a group. The analysis of Kochanek et al. 
(2006), coupled with our findings for the remaining time delay 
systems (77 ~ 2), suggests that strong gravitational lensing time 
delays should be seen as placing valuable constraints on galaxy 
and group formation scenarios, reflecting the interplay between 
baryons and dark matter on these scales. One step to assess the 
impact of galaxy environment on time delays would be to in- 
vestigate the lensing properties of fully simulated groups - are 
the profiles derived from simulated lensing time delays consis- 
tent with these? Such a study would allow us to address the 
impact of complex environments (e.g. for the systems we ex- 
cluded from our sample) on galaxy profiles. 

Ideally, more time delay measurements should be made to 
reduce the errors, and measurements made to increase the accu- 
racy of the astrometry, particularly in systems where the lens- 
ing galaxy's position has notable uncertainties. Moreover, in- 
creasing the sample size of time delay systems is crucial, as it 
puts further constraints on density profiles and reduces the ran- 
dom uncertainties due to large-scale structure. Looking towards 
the future, GAIA will uncover roughly fifty times more multi- 
ply imaged quasars than we have identified today (see Gaia 
study report); of these a few percent will have variable radio 
emission and hence be potential candidates for radio time de- 
lay systems. 
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